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Abstract. The quantum entanglement E of a bipartite quantum Ising chain is com- 
pared with the mutual information / between the two parts after a local measurement 
of the classical spin configuration. As the model is conformally invariant, the entan- 
glement measured in its ground state at the critical point is known to obey a certain 
scaling form. Surprisingly, the mutual information of classical spin configurations is 
found to obey the same scaling form, although with a different prcfactor. Moreover, we 
find that mutual information and the entanglement obey the inequality / < E in the 
ground state as well as in a dynamically evolving situation. This inequality holds for 
general bipartite systems in a pure state and can be proven using similar techniques 
as for Holevo's bound. 
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Figure 1. Quantum spin chain with periodic boundary conditions. A measurement 
of spin orientations yields a classical spin configuration, as indicated by the arrows. 
Dividing the system into two fictitious parts A and B, the respective spin configurations 
will be correlated, expressed in terms of the mutual information 7 A:B . 



1. Introduction 

Recently the study of quantum spin chains from the perspective of quantum information 
theory attracted considerable attention. This applies in particular to entanglement 
studies of quantum spin chains in their ground state p = \ijjo)(' l Po\- In these studies a 
quantum chain is fictitiously divided into two parts A and B (see Fig. [TJ). The quantum 
entanglement between the two parts is then given by 

E A:B = S[p A ] = S[ PB ], (1) 

where S[pa,b] = — tr [pa,b In Pa,b] denotes the von-Neumann entropy of the reduced 
density matrices pa = tr#[p] and p B = tT A [p\- The entanglement E A . B is particularly 
interesting to study in the context of critical quantum chains with an underlying 
conformal symmetry p], which are characterized by long-range correlations in the 
ground state. Using methods of conformal field theory it was shown in [2H6] that the 
entanglement in such systems obeys the scaling form 




where L is the total length of the chain, I is the length of section A, and / is a scaling 
function. In this expression the constant a is non-universal, meaning that it depends 
on the specific microscopic realization of the respective model, while the constant c 
turns out to be universal. Remarkably, c is equal to the so-called central charge of the 
underlying conformal field theory which labels the universality class. For example, the 
Ising universality class is characterized by the central charge c = 1/2. 

To determine the entanglement experimentally one has to perform a variety 
of repeated measurements acting simultaneously on all spins in one of the sectors. 
Such as task is usually difficult to perform. In fact, it would be much simpler to 
measure individual spin orientations locally and to study the resulting classical spin 
configurations. Such a measurement is expected to destroy the existing entanglement 
and to convert it to some extent into classical correlations between the two parts of the 
system. Such classical correlations are usually quantified by the mutual information 

^a : b = H A + H B — H A b i (3) 

where the terms on the right hand side denote the Shannon entropy of the classical spin 
configurations after the measurement in the sections A, B and in the whole chain. 
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Figure 2. Entanglement E A:B (L,l) (squares) and mutual information I A:B (L,l) 
(circles) in the ground state of the quantum Ising model with L = 24 sites and periodic 
boundary condition at the critical point g = 1. The solid and dashed line show the 
scaling functions according to the scaling forms (|9]) and (jlpp . respectively. 

In this paper we address the question to what extent the quantum entanglement 
and the classical mutual information are related to each other. The main results are: 

• In the quantum Ising model the mutual information obeys the same scaling law (J2) 
as the entanglement at the critical point, although with a different prefactor. 

• The mutual information obeys the inequality I A[B < E A:B for general bipartite 
systems in a pure state. 

The paper is organized as follows. In the following Section we first discuss the example 
of the quantum Ising chain, determining the entanglement and the mutual information 
in finite-size systems by different methods. In Sect. 3 we give a general proof of the 
inequality J A:B < E A , B . Finally, in Sect. 4 we present a summary and discuss about the 
relation between our inequality evoked by local measurements and the monotonicity of 
the quantum relative entropy. 

2. Entanglement and mutual information in the quantum Ising chain 

The quantum Ising chain with the length L is defined by the Hamiltonian 

L L 

H = -5>^h -a^of, (4) 
i=i i=i 

where a x ' y,z are Pauli matrices, g is the strength of the transverse field. We use periodic 
boundary conditions by setting L + 1 = 1. The ground state of this model exhibits an 
order-disorder phase transition at the critical point g c = 1: At g = all spins are aligned 
in one direction of the z axis, and as g is increased the magnetization is weakened and 
vanishes at g — 1, which means that the system goes to the paramagnetic phase in the 
z direction at g > 1. 

Scaling behavior E A:B and 7 A:B in the ground state: 

To obtain the mutual information of classical spin configurations in the ground state, 
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we have to compute the probability for each classical configuration after a measurement 
of local spin orientations. If we choose the z axis for the spin orientations, measuring 
erf , erf, . . . , cr£, these probabilities can be computed from the coefficients of the ground 
state wave function |-0 O ) represented in the product basis of the eigenvectors of of. For 
example, if we consider an Ising chain with two spins in the ground state 

lV>o) = «i| tt) + a 2 \ ||) + <z 3 | It) + a 4 \ ||), (5) 

where | "||) = | |)i ® I I) 2 is the product basis of the eigenvectors | f), and | |), 
of o?, the probability P({a z }) to obtain the classical configuration {cr 2 } is given by 
P(tt) = «? , -P(tl) = «2 > -P(lt) = and P(||) = a|. To determine the ground state 
of the Hamiltonian (T4J) at finite L, we use the Lanczos method [7]. 

Each classical configuration obtained from the measurement is now divided into 
two segments, {o~ z } = {o~%,o~ B }. Determining the Shannon entropy of the probability 
distribution in the segments and in the whole chain, 

H A = -J2P(i<})^P({<}), (6) 

{° a) 

Hab=~ Yl P(W Z A^B})^P({a^a B }), (8) 

where P{{a\}) = £ {oS} P(M,o£}) and P({a B }) = E{ CTi } one can 

obtain the classical mutual information J A:B (L, /) = H A + Hb — Hab for a given length / 
of section A. 

On the other hand, to compute the entanglement E A:B (L,l) between sections 
A and B, we perform the partial trace on the density matrix p = l^oX^ol obtained 
from the Lanczos method, leading to the reduced density matrices pa and ps- Since the 
entanglement of a pure state is given by the von- Neumann entropy of the reduced density 
matrices (see Eq. (CQ)), we can determine the entanglement E A . B (L,l) by numerically 
diagonalizing pA or p B - 

Using this method, we calculate E A[B (L, I) and J A:B (L, /) for various values of g, up 

to L = 24. As expected, at the critical point g = 1 the entanglement E A . B (L, I) is found 

to obey the scaling form ([2j) as shown in Fig. [2] with /(£) = - sin(7r£). Surprisingly, the 

mutual information J A:B (L, /) also obeys the same type of scaling form. To illustrate this 

finding, we plot two scaling forms 

1 , / L lir\ , , 

E A:B (L,l =a+-ln -sin— , (9) 

O V 7T L J 

I A:B (L,l) =a' + ^\n(-sm l -pl, (10) 

as a solid and a dashed line in Fig. [21 together with the numerical data obtained by the 
Lanczos method. Here, we have used a ~ 0.478, a' w 0.329, and d ~ 0.715. This means 
that the initial entanglement and the mutual information after the measurement differ 
only by a factor at the critical point. 
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Figure 3. Entanglement E A:B (L,l) (solid line) and mutual information I A:B (L,l) 
(dashed line) in a time-dependent pure state of the quantum Ising model with L = 10 
and I = 5. Initially all spins are aligned in the up direction in z axis. Results for 
g = 0.1 are shown in (a) and for g = 1 in (b). 



Note that the classical mutual information after the measurement is always smaller 
than or equal to the initial entanglement in the ground state. The same observation 
holds in the off-critical case g ^ 1 (not shown here). 

Behavior in a non-stationary pure state: 

The inequality 7 A:B < E A[B is valid not only in the ground state but also in time- 
dependent pure states of the quantum Ising chain. The time evolution of the density 
matrix is given by 

p(t) = e-* m mO))(m\e lHt , (11) 

where H is the quantum Ising Hamiltonian (jlj) and |^(0)) is the initial state. Since p(t) 
is also a pure state, one can obtain E A:B (L, I) and I A[B (L, I) using the same methods as 
in the ground state. In Fig. [3]we show results for g = 0.1 and 1.0 in a chain with L — 10 
and 1 — 5, using an initial state that is fully magnetized in z-direction. As can be seen, 
the quantities oscillate irregularly but satisfy the inequality I A:B < E A[B at any time. 

It turns out that this inequality holds generally for arbitrary local measurements 
on entangled pure states in bipartite system, as will be proved in the following section. 

3. Proof of the inequality J A . B < -E^.b 

To prove the inequality between mutual information and entanglement we use similar 
techniques as for the proof of Holevo's inequality [9] . This suggests that both inequalities 
may be closely related or even equivalent. 
Measurement : 

In the following let us consider an arbitrary bipartite system on T-L AB = 'H A ®'H B which is 
in a pure state p AB = \tjj){tjj\ so that the von- Neumann entropy 5[/> A b] = — tr [p AB lnp AB ] 
vanishes. In such a pure state the entanglement between the subsystems A and B is 
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given by 

E A[B = S[p A ] = S[p B ] , (12) 

where p A = tr B p AB and p B = tr A p AB denote the reduced density matrices. On both sides 
let us now perform local projective measurements 

M A = Y,a\4>a)(^ a \ = $>n a , M B = J2b\4>b)(fo\ = J2 bUb ( 13 ) 

a a b b 

with the projectors Il a = \<f> a ){<f>a\ and lib = |0b)(0b|, which may be thought of as spin 
configuration measurements as in the example given above. This measurement converts 
the pure state p AB into a mixed state 

Pab = ^abpAB^-ab , (14) 

a.b 

where H a b = n a ®rib. Moreover, the measurement completely destroys the entanglement 
between the two subsystems. To see this, let us consider the entanglement of formation 



p' AB = J2q k \k)(k\\ (15) 



E F (p AB ) = Mr£q k S[tr B \k){k\ 

defined as the infimum of the entropy over all possible statistical ensembles represented 
by p' AB . Since E F (p' AB ) is non- negative and the entropy of the particular representative 

P ' AB = ^n a6 tr[n a6P 'ru] (16) 

ab 

vanishes, we can conclude that Ef(p' ab ) = after the measurement. This means that 
the measurement destroys the original entanglement and converts it to some extent into 
classical correlations which can be quantified by the mutual information J A:B in Eq. ([3]). 
Note that the classical mutual information J A . B is equivalent to the quantum mutual 
information I(p' AB ) of the post-measurement state p' AB , i.e., 

Ja,b = I{p AB ) = S[p' A ] + S[p' B ] - S[p' AB ] , (17) 
where p' A = tr B [p AB ] and p' B = tr A [p AB ]. 



Expressing the measurement as an isometry in an extended space: 

According to the Stinespring theorem [5] the measurement process fll6p can be carried 
out by embedding p AB in a higher- dimensional Hilbert space, performing a unitary 
transformation on it, and finally tracing out the additional degrees of freedom. To this 
end let us extend the Hilbert space % AB = 'H A <S)'H B by an auxiliary space "H AB = 'H A ®'H B 
whose task will be to store the measurement outcome encoded in the form of orthonormal 
basis vectors \<p a b) = 1 0a) <8> | </>&)• Moreover, let us define the linear map 

v ■u AB ^n AB ®u AB : \^)->v\4j) = J2 

ab 
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which obeys the condition VV^ = 1 so that it maps the initial state isometrically onto a 
subset of the extended space like a unitary (entropy-preserving) transformation. With 
the corresponding extended density matrix 

Pabab := Vp AB V^ = y^y^(0afe|PA B |</W) (\4>ab)(<Pa'b'\<^\4>ab)((f>a'b'\) (19) 

ab a'b' 

the measurement process (Tl6|) can now be written as 

p' AB = tr AB [Pabab] • (20) 
Schmidt decomposition: 

According to Schmidt's theorem, any quantum state G % AB can be decomposed into 

= ^\/A^a>®Kb> (21) 

i 

with certain vectors |z A ) G "H A , \i B ) G H B and probabilities Aj G [0, 1], called Schmidt 
coefficients, which sum up to 1. This means that the initial state can be written as 

p AB = |-0> <VI = $^\/Vv Na)(Ja| ® Kb)0'b|- (22) 

Calculating the partial traces we obtain 

p A = tr B [p AB ] = ^2 V 7 Na>0a| tr B |z b >0b| = X^'^ A ^ A I ( 23 ) 

=<5;j 

and a similar expression for p B , meaning that the Shannon entropy of the Schmidt 
coefficients is equal to the initial entanglement: 

E A . B = -J2 X i lnX i ( 24 ) 

i 

Encoding the Schmidt decomposition in another auxiliary space: 

Let us now introduce another auxiliary Hilbert space "H c with canonical basis vectors 
\ij) whose task will be to store pairs of Schmidt coefficients. On the combined Hilbert 
space l-L AB c = Hab <8> He we define the density operator 

Pabc := I z 'a)0a| ® |«b>0'b| ® (25) 

ij 

If we compare this operator with its own diagonal part 



Pabc 



^ \ |*a)(*a| <S> |« b )(*b| ® (26) 



it is easy to check that integer powers of these operators have always the same trace, 
i.e. tr[p ABC ] = tr[p ABC ] for all k G N, meaning that they have the same entropy 

SVabc] = SIPabc] = - ^2 K In Ai = E A:B ■ (27) 
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We now use the map V denned in Eq. (lisp to extend this operator even further to the 
space Hababc = %ab <S> T-Lab <8> "He by defining 

Pababc = (V ® l c ) Pabc {V ] <g> l c ) 

= XX^X^ {<f>a\iA){3A\<j>a?){<l>b\iB)(jB\<h>) (28) 

ab a'b' ij 

X \(pab)(4>a'b'[ ® \^ab)(<Pa'b'[ ® K^Kuj " 
AB AB C 

As V is an isometry, this extension does not change the entropy, hence 

S[p A BABc] = E A -.B ■ (29) 

Tracing out the original Hilbert space "H AB and the auxiliary space H c '■ 

Tracing out the original Hilbert we have ti AB [\(j) ab )((j) a i b /\] = 5 aa >5 bb >, leading to 

PASO = tr ab [PaBABc] = ^EV^i (<Pa\iA)(jA\<Pa)(<Pb\iB)(jB\<Pb) (30) 

ab ij 

X \<Pab)((pab\ <S> ■ 



Since the density matrix in Eq. (1301) has only diagonal elements, the von-Neumann 
entropy of pabc and its reduced density matrices can be obtained easily as follows: 

SVabc] = - ^2 ^ ln ^ ~ ^^(V^j (0tt|iA)(7A|0a)(06|*B>(7B|06) 
i ab ij 

x \n{(p a \i A ){j A \(f) a ){4) b \i B ){j B \(f) b ) ) , (31) 



S[p AC ] = - Aj In \j - Xj {(j)a\i A ){i A \(j)a) In {(j)a\i A ){i A \^a) , (32) 

i a,i 

S[Pbc] = -^XilnXi -^Xi (<p b \i B )(i B \(p b ) ln((f) b \i B )(i B \(j) b ) , (33) 

i b,i 

where p AC = trg[pABc] an d p B c = tr^ [pabc] • Combining Eqs (13T]) - (133|) . one obtains a 
useful relation 

S[p AC ] + S[p BC ] - S[p ABC ] = E A: b ■ (34) 

Tracing out the auxiliary space "H c we obtain the reduced density matrices p AB ,p A , and 
p B in which the outcome of the measurement is stored. Therefore, the classical mutual 
information after the measurement can be expressed in terms of these density matrices 
by 

S[pa] + S[p B ] - S[p AB ] = I A:B . (35) 
Apply strong subadditivity: 

To prove the inequality J A :b < E A:B , we use the strong subadditivity of the von-Neumann 
entropy [10] which is known to hold in the two following equivalent forms |llj : 

S[p ABC ] + S[Pb] < S[p AB ] + S[p BC ] (36) 

S[p A ] + S[p B ] < S[p AC ] + S[p BC } . (37) 
Using the first relation the mutual information can be written as 

Ia-.b = S[p A ] + S[p B ] - S[p AB ] < S[p A ] + S[pbc] - S[p ABC ] (38) 
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Replacing A -B- B one obtains in the same way 



h:B < S[ps] + S[p Ac \ ~ SVaBc] • 



(39) 



Adding ( 1381) and ( 1391) gives an inequality where we again apply the second version of 
strong subadditivity (I37I) 



proving the initial assertion. 
4. Conclusion 

We have investigated the relation between the initial entanglement and the classical 
mutual information after local projective measurements in a bipartite quantum system 
which is initially in a pure state. For the quantum Ising chain with periodic boundary 
conditions, which is divided into two segments, it is found that both quantities obey 
the same scaling form at the critical point, differing only in the prefactor of the scaling 
function. For the entanglement we find the prefactor c/3 = 1/6 in agreement with 
predictions of conformal field theory, whereas the prefactor of the mutual information 
is found to be less than 1/6. 

Furthermore, we have observed that classical mutual information cannot exceed 
the initial entanglement in arbitrary pure states. The inequality J A:B < E A[B has 
been proved generally by successively expansions and reducing the Hilbert space and 
by applying strong subadditivity of the von-Neumann entropy. For entangled pure 
states, we conclude that a local projective measurement destroy the original quantum 
correlation and converts it into classical mutual information bounded from above by the 
initial entanglement. 

In a general case, it is more convenient using the quantum mutual information 
HPab) = S[Pa] + S[Pb] — S[p AB ] to measure the correlation in a quantum system. 
According to the monotonicity of the quantum relative entropy, the local projective 
measurement does not increase the quantum mutual information [12], that is, 



where p' AB is the post measurement density operator. If the initial system is in a pure 
state, Eq. (jUJ) yields J A:B < 2E A:B which is consistent with our inequality. As a future 
work, generalized inequalities induced by the measurement on mixed state will be studied 
in a context of the quantum mutual information. 
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